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The Dirac equation for the Coulomb problem is restated by incorporating 
a nonlinear effective interaction into the Dirac Hamiltonian: one keeps the 1/r 
dependence for the Coulomb field, but the coupling constant is modified by 
a factor depending on the n {principal quantum number) power of the mean 
value of the Hamiltonian. In this simple context we study the Lamb shift 
and the hyperfine splitting of the s-levels of hydrogenic atoms. We discuss 
to what extent the corresponding calculations fit the energy splittings to the 
appropriate order in the fine structure constant. 

PACS number(s): 03.65.P, 32.10.F, 31.30.J 

I. INTRODUCTORY BACKGROUND 

In 1947, Lamb and Retherford observed that the 2si/2 and the 2pi/2 energy levels 
of the hydrogen atom were split: the 2pi/2 energy level was depressed more than 1000 MHz 
below the 2s 1/2 energy level. The original theory of a Dirac electron in a classical Coulomb 
potential predicted that the energy levels of the hydrogen atom should depend only on the 
principal quantum number n and the total spin j, so these two levels should be degenerate. 

The calculation of the Lamb-shift is rather intricate, because one is dealing with the 
hydrogen atom (Z = 1) as a bound-state problem, and also because we must sum over 
all radiative corrections to the electron interacting with a Coulomb potential that modify 
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the naive \l/'''y4o\l/ vertex. These corrections include the vertex correction, the anomalous 
magnetic moment, the self-energy of the electron, the vacuum polarization graph, and even 
infrared divergences. 

The original nonrelativistic bound-state calculation of Bethe [Q, which ignored many of 
these subtle higher-order corrections, could account for about 1000 MHz of the Lamb-shift, 
but only a fully relativistic quantum treatment could calculate the rest of the difference. 
To begin the discussion, we first see that the vacuum polarization term can be attached to 
the photon line, changing the photon propagator. This, of course, translates into a shift in 
the effective coupling of an electron to the Coulomb potential |Q. We know from ordinary 
nonrelativistic quantum mechanics that, by taking matrix elements of this modified potential 
between Coulomb wave functions, we can calculate the first-order correction to the energy 
levels of the one-electron atom due to the vacuum polarization graph. 

This discussion can be generalized to include the other corrections to the calculation of 
the Lamb-shift. The method is the same: corrections to the vertex function \I^^7^\1/, take the 
zeroth component, and then take the low-energy limit. If we add the various contributions 
to the vertex correction, we find the well-known formula: 

(1) 

to order a (Za)^ In (Za) , where m is the mass of the electron and a = e^/hc is the fine 
structure constant. The term 
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is also known as the Bethe logarithm which has to be evaluated numerically. In the above 
K = ± (j + 1/2) for / = j ± 1/2, and is a nonrelativistic atomic state. 

The vertex correction, for Z = 1 for example, gives us a value of 1010 MHz. The 
anomalous magnetic moment of the electron contributes 68 MHz. The vacuum polarization 
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graph contributes —27.1 MHz. Adding these corrections together, we find, to the lowest 
loop level, the Lamb-shift to within 6 MHz accuracy. 

Since then, higher-order corrections have been calculated, so that the difference between 
experiment and theory has been reduced to 0.01 MHz . Theoretically, the 2si/2 level is 
above the 2pi/2 energy level by 1057.864 MHz. The experimental result is 1057.862 within 
0.02 MHz. This is an excellent indicator of the basic correctness of QED 

Another important effect, not contained in the energy levels of the Dirac-Coulomb prob- 
lem, arises from the interaction between the magnetic moment of the nucleus and the mag- 
netic moment of the electron. In the case of the hydrogen atom, for instance, when we 
combine the electron spin with the proton spin, the net result is S" = 1 (triplet) or S" = 
(singlet), where S is the quantum number corresponding to the total spin. Since the magnetic 
interaction is dependent on the relative orientation of the two magnetic dipole moments, each 
level of the hydrogen atom characterized by njl, is further split into two sublevels corre- 
sponding to the two possible values of S even in the absence of any external magnetic field. 
This is known as the hyperfine splitting [pHTlJ . 



Using the nonrelativistic wave function \l/nK^'', we obtain the energy shift 

^E^^''^ = ^g, ^ i6s,i - 36s,o) 5.,-,a\ (3) 

where gp = 2 (1 + Hp) = 5.58568 is the (yf-factor of the proton. Note that the order of 
magnitude of this splitting is the fine- structure splitting multiplied by m/Mp. For the 2si/2- 
state, the above energy difference corresponds to a radio microwave of 1420 MHz. It is an 
accurately measured quantity: 1420.40575180 within at least 3 x 10^^ MHz |jl3|. 

There are other corrections to the Dirac formula. First, we must take into account the 
motion of the nucleus since the mass of the nucleus is not infinite. A major part of this 
correction can be taken care of if we use the reduced mass in place of m. Second, the finite 
size of the nucleus, especially for the s states which are sensitive to small deviations from 
Coulomb's law at close distances: in the interesting case of the 2s state of the hydrogen 
atom, however, we can estimate the energy shift due to this effect to be only 0.1 MHz. 

The utility of Dirac theory in atomic physics is not limited to light hydrogen-like atoms. 
For heavy atoms where Za is not very small compared with unity, the relativistic effects 
must be taken into account even for understanding the qualitative features of the energy 



levels. Although we cannot, in practice, study one-electron ions of heavy atoms, it is actually 
possible to check the quantitative predictions of Dirac theory by looking at the energy levels 
of the innermost electrons of high Z atoms which can be inferred experimentally. Similar 
studies have been carried out with muonic atoms [|1^ . 



This paper concerns a simple, and restricted, approach to the study of the splitting 
structure of the Dirac-Coulomb energy levels, given in Eqs.(|l]) and (|]), in the context of 
relativistic quantum mechanics. This procedure circumvents second quantization on both 
the electromagnetic and the electron fields. However, it does not pretend to be an alternative 
way of reproducing the corresponding detailed and precise modern calculations by QED. In 
fact, for the case of the Lamb-shift, we do in fact make use of second quantization information 
(the self-energy of the electron, the vertex correction, the anomalous magnetic moment, and 
the vacuum polarization graphs to order a iZa)^ ni(?^ to define the basic coupling constant. 
In Sec. 11, we present the general approach. It considers the introduction of a nonlinear 
effective interaction into the Dirac Hamiltonian. Although the radiative processes involved 
in the Lamb-shift will somehow be hidden in the corresponding effective Hamiltonian, it is 
illustrative to view them from a different perspective. The case of the hyperfine structure is 
treated in the same framework, though it does not involve radiative corrections to the first 
level, i.e., when only relativistic corrections are taken into account. 

Higher order corrections to the Lamb-shift to order a{Za)^mc'^ and to the hyperfine- 
splitting are considered in Sec. HI. Finally, Sec. IV contains our conclusions and some open 
questions. 



II. EFFECTIVE NONLINEAR COULOMB-LIKE INTERACTION 

In this paper we want to show that the splittings (|ip and can be derived from the 
Dirac-Coulomb problem by incorporating a nonlinear effective interaction into the Dirac 
Hamiltonian. To this end, we shall assume that the interaction does not modify the ^ (l/o') 
{q = |q|) law in the Coulomb-like gauge (central) field. This radial behavior is one of the 
few cases where the Dirac wave equation can be solved analytically. The problem to study 
is then the following 
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in the standard Dirac representation 113]. Here the normahzation condition is 



and 



Af^ (g) = % 



Z\e\ 



(4) 



(5) 



(6) 



is an effective potential, where is a diagonal constant matrix: It is taken not to be 
g-dependent since we shall not consider corrections in the vicinity of the nucleus. The 
dimensionless elements of this matrix depend on the quantum numbers represented by 
labeling the stationary states \l/£;^(q, t) = exp (— i_E't,t/fi) \I'£;^(q), whose concrete structure 
will be specified, in each case, below. 

Before solving the eigenvalue problem associated with (H), we recall that the operators 



with L= qxp the orbital angular momentum operator, are constants of motion: 
0, 



(7) 



H,K 



0. Following a standard procedure [15], the stationary states of energy E can 



be written as 



^E.(q,t) 



^a(q,t) 



where yjj^i are the normalized total angular momentum functions, with 



K y. 



hny. 



(8) 



(9) 



where la = j i: 1/2, = j =p 1/2 when k = ± (j + 1/2) . Let us write g^, in the general form 

(10) 



lu = lguAi + P) + lguAi-P)- 



2 " ' ' ' 2 

where Ua^h are labels defined in Eq.(^). Thus the Dirac equation is equivalent to the set of 
first-order (nonlinear) differential equations 

ca-pV^b(q) = (^E^ - mc^ + S'^a^^ V'a(q), ccr-p^/^a(q) = + mc^ + g'^s^^ V^b(q)- 

(11) 



From d), d) and (|TTD we find that 



->('■)) (12, 

where RA,B{r) = ripa,b{r) and 

r = A/^i^2g, Ml = — , M2 = — , Za = Z— . (13) 

Next we look for solutions in the form of series 

-RaI?") = exp(-r)r* ^a^r^, Rsir) = exp{-r)r'''^b^r'^. (14) 
^J.=o /x=o 

Thus from ([T^p and ([T^) we get 

(s - K)ao- gu.Zaho = 0, 
^ g^^Zaao + {s + k) bo = 0, 

and 

(s + /i + k) - a^_i - gy^Zahf, - ^/WJMih^^i = 0, 
(s + - k) 6^ - + gz., ^aa^ - a/Mz/Mio^.i = 0, 

Given that ctO) ^0 7^ 0, from (|r^) we obtain 



for /i = 0, (15) 



for n>0. (16) 



^ = ±7 i'^' - iZafs^.Su,) ^ ±sj [k^ - (Zaf) > -\ • (17) 

The negative sign must be excluded since it would make the functions Ra,b singular at the 
origin. Choosing /i = n' + 1 and a„/+i = 6„'+i = 0, in order to terminate the series, we have 



that hn' = —a^i \/~M^2TM\- Then from (|T^) we get 



2 {s + n') VM1M2 = Za {g^M^ - g^^M^) . (18) 
Finally, from (p!7| ) and (|18|) we obtain the energy eigenvalues from 



2 {s + n - \J imc^) - E^^ = Zag^^ [mc^ + En^,) - Zag^^ {mc^ - E^,^) , (19) 
where 



n = n' +\k\ =n' + j + ^ (20) 

is the principal quantum number. Note that for the point nucleus there exist bound solutions 
(for K = -1) only up to Z ^ ^ 1/a. 

Given the fact that < 1/q >,j,(nr)^ ^/n^, we expect that {g^^ — 1) a 1/n, to be able to 
meet the factor in both (^ and (|^). Thus, in the examples that follow, we shall choose 
a particular form for g^,^ and g^,^: 



= 1 - Xn^i^a, Za) (^1 - Hd >^(^^^--)) ) , (21) 
nKib = 1 - XnKh{(^, Za) ^1 - Hd >q,(oi^ 



bVa — gnnla 

gub = griKii, = J- - '^nKl^,\U^, ZjtX) ^± - ud ^^(Dirac) 

where Hd = Hd/ttlc^. The factors XnKi{(^, Za) may depend on the original Dirac-Coulomb 
potential quantum numbers Ua^b = nula^b , the binding (powers of Za), and radiative cor- 
rections (powers of a). The term containing the expectation value of Hd has the required 
property, namely 

(l - (< Hd >^(-^c))'^) = H (Za)' + O {{Za)') a i. (22) 

Thus the factors gnni includes a nonrelativistic coupling and a radiative term: XnKi{c(,Za) 
(see Eq.(pl|)), times a purely quantum relativistic factor: ^1 — ^< Hd >^„(Dirac) 
Now from ([T9|) and (^TJ) we finally get 



2{s + n- \k\) J{mc^)' - E^, = 2En.Zagnni^ + O {mc^ (Za)' (g„«i^ , - l)) . (23) 



On the right-hand side of (|23D we are neglecting a term proportional to (Za)^ (gnK«„ — l) 
which is, for instance, of order a (Za)^ for the case of the Lamb-shift. 



A. The Lamb-shift 



As a first instance, we shall consider an approximation to the Lamb-shift. To this 
end notice that < Hd >^^{Dirac)= Ennig = I), with Hd = H{g = I) the Dirac-Coulomb 

(Dirac) 

Hamiltonian, where the eigenstates satisfy the normalization condition 

/ <r"(q)<r"(q) A=l. (24) 
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Let 

A„«K«, Za) = ^a(^ [l^i + ^ - 2 In {Za)^ 5i,o + l ^^2\l\-l) " ' ^^^^ 

The effective gauge potential to be considered is of the form 

Af\q)=t^r'^^, (26) 

where 

Sitr^ = gitrH«, Za) = l- A„.,(«, Za) (l - (< Hn >^(-a.))") . (27) 

Thus g^nK™^^ diminishes the Coulomb binding —Za/q, and as a consequence, the s levels are 
pushed higher. From (p!9| ) and we find 

is + n-\K\) {mc-f - [Et-'^y = Zag^^'^E^r'^ + O {mc'a (Zaf) . (28) 

Expanding this equation in powers of a, we find the spectrum 

(4-' - - iZ. f H- (2i - li^) iZa)' (29) 



We observe that the last term in (|29|) corresponds to the Lamb-shift splitting to the Dirac 
levels to order a{Za)^ p|. The Coulomb potential for the hydrogen atom {Z = 1) and the 
corresponding Coulomb-like radiative corrections for n = 4, 8, 12, are shown in fig. L 



B. The hyperfine splitting 

The second instance regards the hyperfine splitting in the energy levels of the hydrogen 
atom {Z = 1). In hydrogenic atoms, the interaction of the magnetic moment of the orbital 
electron with the magnetic moment of the nucleus leads to a splitting of the fine structure 
levels with fixed orbital angular momentum / (k = — 1) into the hyperfine structure levels. 
In this case we choose 

gi«r^ = 1 - A., (l - (< if^ >^,,^^^^, ) ") , (30) 
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where 



^ki = -gp 



m 



{6s,i - SSsfi) 5z,o 



(31) 



with Qp = 2 [1 + Kp) , Hp = 1.79284, the (/-factor of the neutron. Thus from dTPj ) and 

we get 



2{s + n 

from which we get 

( /?(Hyp) _ 



(Hyp)V 



FA ■ 



(32) 



mc ) mc = ^ + 
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1 1 





1 1 
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f m\ 


(^5,1 - 




(Mp) 


35 



a 



(33) 



+0 (a^) . 

This is again the correct spectrum for the relativistic Dirac levels to order O (a'^) |1^ . 

Notice that the first term in both ( p9D and ( p3D gives the energy spectrum of the bound 
states in the non-relativistic approximation. The second terms are the corresponding lead- 
ing corrections to the Balmer formula: the fine splitting. These expressions are a conse- 
quence of the modification Za gnuZa introduced by the nonlinear factor ^< H >^(Dirac)^ 
(together with A„k(«5 Za) in the case of the Lamb-shift), which depend on the solution of 
the original {gnu = I) eigenvalue problem itself. 

Finally, an interesting exercise to consider regards the eigenvalue problem containing 
both corrections: the Lamb-shift and the hyperfine splitting for the case of the hydrogen 
atom {Z = 1). This is easily done by defining 



;^(Lamb) 



+ 



?(Hyp) 



(34) 



Replacing g^^ in (^) and expanding in powers of a yields 

/i^ ■2\ I 2 ICt^ /'SI 111\4 1 f m\ 1 oe-\e- 

- mc ) /mc = --_ + (^-- - --_ j a + -j, j -j (S,, - SSs.o) K-.a (35) 



k{2\k\ - r 



(1 - S^-i) a 



+0 (a^) , 

which is the spectrum of the hydrogen atom to order a^. Notice that any series expansion 
of In (a) rapidly overcomes the relativistic (Dirac) corrections of Enn beginning from order 
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III. HIGHER ORDER CORRECTIONS 



To achieve a still better quantitative agreement, several contributions of higher order 
must be included namely corrections to order mc^a [Zaf , mc^a [Zaf and m(?cP' [Zaf . 
Here we only want to examine the mc^a [Zoif term (the 2nd order binding) as the next 
correction to the Lamb-shift: 2si/2 — ( p9D for the hydrogen atom (see fig. 2). To 

this end we can approximately fit the corresponding QED theoretical value by making the 
replacement 

XnKi{a, Za) -> (1 + Za) Xn^iia, Za), (36) 
in the A„k« factors contained in (127|) . For Z = 1, the correction calculated with (|T9|) , (|27D 



and (|36|), and the QED value 0, both to order a {Zaf , are respectively: 

bE ^ 7.663 MHz , (5^)(qed) - 7-243 MHz , (37) 

which corresponds to a ^ 5% of discrepancy. These calculations have been performed with 
the value a ^ 1/137.036 ^ 7. 29735 x 10~^ for the fine structure constant. 

On the other hand, in fig. 2 we observe that there is good agreement between the 
present approach and the QED numerical values for Z = 1,2, 10. In fact, the Lamb-shift 
for Z = 1, n = 2, is given by 



^^^Lamb) ^ 04Q_54 MHz , (^AE^^^^""^^ = 1, 046.45 MHz , (38) 

where 

^^(Lamb) ^ ^(Lamb) (^^gy^) - E^^^'^^\2py2) , (39) 

(see the red line in fig. 2); while for Z = 10, n = 2, we get 

^^(Lamb) ^ g2 X 10^ MHz , ( AE^^'^^'^A = 4, 860.51 x 10^ MHz , (40) 

V / QED 



with a ^ 8 % of discrepancy. 

For the hyperfine structure we make 



Air^ - X^S'"^ («, Za) = (1 + Veit) + S {a, Za)) A^f , (41) 
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where (5(Breit) = (3/2) (Za)^, and S{a,Za) contains higher order corrections: radiative, 
binding, finite mass and structure of the proton Higher (radiative and binding) 

order corrections to the Lamb-shift seem to be meaningless and difficult to reproduce by 
using this simple approach . 



IV. OUTLOOK 



In this article we have introduced a nonlinear effective Coulomb-like gauge potential 
(electromagnetic in nature) into the Hamiltonian of a Dirac particle to describe the Lamb- 
shift and the hyperfine structure of the energy spectrum. In all the calculations we have kept 
track of the corresponding orders in both the fine structure constant a and the combination 
Za introduced by the g^^ matrix factor. In Eq. (p!9| ) there still is room for further research. 
For instance, we can study systems (apart from the example of the Lamb-shift) where gu{ct) 
is of order ^ > (^ some integer number) in a and Za. In addition to this, the prospect 
of defining an iterative recursion procedure from (|19[) and (^) in order to improve the 
(eventually convergent) coupling constant Za should be studied. 

Much work needs to be done to give a precise characterization of the various physical 
models included in Eq.(|l^). Particularly, there still is lacking a thorough explanation of the 
structure of the gy(a;) factors contained in the different effective Hamiltonians that we have 
considered. 

This work was supported by Direccion de Investigacion, Universidad de Concepcion, 
through grant #96.011.019-1.0, and Fondecyt through grant #1970995. 
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Figure captions 



Fig. 1. Schematic plots of the Coulomb potential per unit charge y4o(g)/|e| for the 
hydrogen atom and the corresponding Coulomb-like (radiative) |A„K/^(a, Za)\ Ao{q)/ \e\ {k = 
— 1, la = 0) corrections (Re) for n = 4,8, 12. 

Fig. 2. The lamb-shift 2si/2 — for Z = 1, ...,40. The calculations with the present 
approach to order a (Za)^ (green color) are improved with the corrections to order a [Zaf 
(red color). The continuous lines are drawn in order to visualize the corresponding patterns. 
There is good agreement for values of Z in the range 1, .., 10 with those of QED (see Refs. 



Tql ) with a relative accuracy between 0.01 %. {Z = 1) and 8 % {Z = 10). 
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